Presented at the workshop: 
"Fundamental Physics of Ferroelectrics" , 
Aspen Center for Physics, February 13-20, 2000 



What makes an insulator different 
from a metal? 

Raffaele Resta 



INFM - Dipartimento di Fisica Teorica, 
Universita di Trieste, Strada Costiera 11, 1-34014 Trieste, Italy 



Abstract. The insulating state of matter is characterized by the excitation spectrum, 
but also by qualitative features of the electronic ground state. The insulating ground 
wavefunction in fact: (i) displays vanishing dc conductivity; (ii) sustains macroscopic 
polarization; and (iii) is localized. The idea that the insulating state of matter is a 
consequence of electron localization was first proposed in 1964 by W. Kohn. I discuss 
here a novel definition of electron localization, rather different from Kohn's, and deeply 
rooted in the modern theory of polarization. In fact the present approach links the two 
features (ii) and (iii) above, by means of essentially the same formalism. In the special 
case of an uncorrelated crystalline solid, the localization of the many-body insulating 
wavefunction is measured — according to our definition — by the spread of the Wannier 
orbitals; this spread diverges in the metallic limit. In the correlated case, the novel 
approach to localization is demonstrated by means of a two-band Hubbard model in 
one dimension, undergoing a transition from band insulator to Mott insulator. 



INTRODUCTION 

The present contribution, dealing with such a general issue as the difference 
between insulators and metals, may appear out of place at a workshop focussing on 
the "Fundamental Physics of Ferroelectrics" . But indeed, the present results are 
in a sense the ultimate developments of the theory of polarization [1-5] based on a 
Berry phase [6], which is crucial to the modern understanding of ferroelectrics [7]. 
The link between the present subject and polarization is simply stated: insulators 
sustain nontrivial bulk polarization, metals do not. 

To the present purposes, materials are conveniently divided in only two classes: 
insulators and metals. We use here the terms in a loose sense as synonym of 
nonconducting and conducting: an insulator is distinguished from a metal by its 
vanishing conductivity at low temperature and low frequency. This qualitative 
difference in the dc conductivity must reflect a qualitative difference in the orga- 
nization of the electrons in their ground state. So the question we are going to 
address is: is it possible to find a pure ground state property which discriminates 
between an insulator and a metal? Before proceeding to answer, let me discuss the 



alternative phenomenological characterization of the insulating/metallic behavior: 
instead of making direct reference to the dc conductivity, we address macroscopic 
polarization. 

Suppose we expose a finite macroscopic sample to an electric field, say inserting it 
in a charged capacitor. Then the induced macroscopic polarization is qualitatively 
different in metals and insulators. In the former materials polarization is trivial: 
universal, material-independent, due to surface phenomena only (screening by free 
carriers). Therefore polarization in metals is not a bulk phenomenon. The opposite 
is true for insulators: macroscopic polarization is a nontrivial, material-dependent, 
bulk phenomenon. We can therefore phenomenologically characterize an insulator, 
in very general terms, as a material whose ground wavefunction sustains a bulk 
macroscopic polarization whenever the electronic Hamiltonian is non centrosym- 
metric. From this definition it is clear that the modern theory of polarization, 
based on a Berry's phase, can lead to a better understanding of the insulating state 
of matter. 

This paper is organized as follows. First we briefly outline Kohn's theory of 
the insulating state [8]. Then we address dipole and localization for a lone elec- 
tron in a Born-von-Karman periodic box; subsequently, we apply similar ideas to 
an extended system of N electrons in a periodic box, addressing macroscopic po- 
larization and electronic localization in the many-body case. We then show how 
this works for a crystalline system of independent electrons, and finally — following 
Ref. [9] — we demonstrate localization for a model correlated system displaying two 
different insulating phases. For the sake of simplicity, we explicitate here the rele- 
vant algebra only for the case of one-dimensional electrons. The generalization to 
three-dimensional electrons can be found in Refs. [10-12]. 

THE INSULATING STATE AND KOHN'S THEORY 

Within any classical theory, the electronic responses of insulators and metals are 
qualitatively described by "bound" and "free" charges, respectively. Microscopic 
models for such charges are provided by the Lorentz theory (insulators) and by 
the Drude theory (metals). Within the former, each electron is tied (by an har- 
monic force) to a particular center; within the latter, electrons roam freely over 
macroscopic distances, hindered only by atomic scattering potentials. Therefore, 
from a purely classical model viewpoint, one explains the insulating/metallic be- 
havior of a material by means of the localized/delocalized character of the electron 
distribution. 

Switching to quantum mechanics, this clearcut character of the electron distribu- 
tion is apparently lost. Textbooks typically explain the insulating/metallic behavior 
by means of band structure theory, focussing on the position of the Fermi level of 
the given material: either in an band gap (insulators), or across a band (metals). 
This picture is obviously correct, but very limited and somewhat misleading. 

First of all, the band picture applies only to a crystalline material of independent 



electrons: a very limited class of insulators indeed. Noncrystalline insulators do 
in fact exist, and the electron-electron interaction is a fact of nature: in some 
materials the insulating behavior is dominated by disorder (Anderson insulators), 
in some other materials the insulating behavior is dominated by electron correlation 
(Mott insulators). Further classes of insulators are also known, as e.g. excitonic 
insulators. Therefore, for a large number of insulators, the band picture is totally 
inadequate. 

Second: even for a material where a band-structure description is adequate, the 
simple explanation of the insulating/ metallic behavior focusses on the spectrum of 
the system, hence on the nature of the low lying electronic excitations. Instead, 
the qualitative difference in the dc conductivity at low temperature must reflect 
a qualitative difference in the organization of the electrons in their ground state. 
Such a difference is not evident in a band-structure picture: the occupied states 
are of the Bloch form both in insulators and in metals, and qualitatively rather 
similar (in particular those of simple metals and of simple semiconductors). 

In a milestone paper published in 1964, Kohn was able to define the insulating 
state of matter in a way which in a sense is close to the classical picture. In fact 
he gave evidence that electron localization is the main feature determining the 
insulating behavior of a many-electron wavefunction [8], thus restoring the same 
basic distinction as in the classical picture: the key is how to define and to measure 
the degree of electronic localization, visualizing in a qualitative and quantitative 
way the peculiar organization of the electrons which is responsible for the insulating 
state of matter. 

As previously stated, a superficial look indicates that electrons are roughly speak- 
ing equally delocalized in insulators and in metals. One needs therefore a sharp 
criterion which singles out the relevant character of the wavefunction. Kohn's cri- 
terion is the following [8,13,12]: the many-electron wavefunction is localized if it 
breaks up into a sum of functions ^ = Yl j which are localized in essentially 
disconnected regions TZj of the configuration space. Any two such \l/j's have ex- 
ponentially small overlap. Under such a localization hypothesis, Kohn proves that 
the dc conductivity vanishes. 



The dipole moment of any finite iV-electron system in its ground state is a 
simple and well defined quantity. Given the many-body wavefunction \1/ and the 
corresponding single-particle density n(r) the electronic contribution to the dipole 
is: 



DIPOLE AND LOCALIZATION FOR A SINGLE 

ELECTRON 
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FIGURE 1. The distribution 1^(^)1 2 of a single-particle orbital within periodic 
Born-von-Karman boundary conditions 

where R = Y2iLi r «- This looks very trivial, but we are exploiting here an essential 
fact: the ground wavefunction of any bound iV-electron system is square-integrable 
and vanishes exponentially at infinity. Going at the very essence, we simplify 
matter at most and we consider in the present Section only a single electron in one 
dimension. The dipole (or equivalently the center) of the electronic distribution is 
then: 



where again we understand ip(x) as a square-integrable function over R. This is not 
the way condensed matter theory works. Because of several good reasons, in either 
crystalline or disordered systems it is almost mandatory to assume BvK boundary 
conditions: the wavefunction ip(x) is periodic over a period L, large with respect 
to atomic dimensions. 

Adopting a given choice for the boundary conditions is tantamount to defining the 
Hilbert space where our solutions of Schrodinger's equation live. By definition, an 
operator maps any vector of the given Hilbert space into another vector belonging to 
the same space: the multiplicative position operator x is therefore not a legitimate 
operator when BvK are adopted for the state vectors, while any periodic function of 
x is legitimate: this is the case e.g. of the nuclear potential acting on the electrons. 

Suppose we have an electron distribution such as the one in Fig. 1. The main 
issue is then: how do we define the center of the distribution? Intuitively, the 
distribution appears to have a "center", which however is defined only modulo 
the replica periodicity, and furthermore cannot be evaluated simply as in Eq. (2), 
precisely because of BvK. 

Solutions to this and similar problems have been attempted several times: many 
incorrect papers — which will not be identified here — have been published over the 
years. The good solution has been found by Selloni et al. in 1987 by means of a 
very elegant and far-reaching formula [14]. According to them, the key quantity 




(2) 



for dealing with the position operator within BvK is the dimensionless complex 
number 3, defined as: 



3 = <#^ X M= [ L dxe^mx)\ 2 , (3) 
Jo 

whose modulus is no larger than 1. The most general electron density, such as 
the one depicted in Fig. 1, can always be written as a superposition of a function 
n\ oc (x), normalized over (—00, 00), and of its periodic replicas: 

00 

|-0(a?)j 2 = ^ n ioc (x-x -mL). (4) 

m=— 00 

Both x and n\ oc (x) have a large arbitrariness: we restrict it a little bit by imposing 
that Xq is the center of the distribution, in the sense that dxxri\ oc (x) = 0. 
Using Eq. (4), 3 can be expressed in terms of the Fourier transform of ri\ oc as: 

3 = e^°n loc (-^). (5) 

If the electron is localized in a region of space much smaller than L, its Fourier 
transform is smooth over reciprocal distances of the order of L~ l and can be ex- 
panded as: 

2tt x 1 /2tt x 2 



nioc(--^) = 1 -2 [t) J dxx n ^ x ) + °( L )■ ( 6 ) 

A very natural definition of the center of a localized periodic distribution |-^(a;)| 2 
is therefore provided by the phase of 3 as: 

(x) = log 3, (7) 

which is in fact the formula first proposed by Selloni et al. [14]. The expectation 
value (x) is defined modulo L, as expected since ^(x)) 2 is BvK periodic. It is also 
worth to observe that for an extremely delocalized state we have ^(a;)! 2 = 1/L and 
3 = 0: hence the center of the distribution (x), according to Eq. (7), is ill-defined, 
as one would indeed expect. 

So far, we have not specified which Hamiltonian we were addressing when dis- 
cussing electron distributions ^(a;)] 2 of the kind depicted in Fig. 1. It is however 
obvious to imagine that the wavefunction ip(x) is the eigenstate of a (periodi- 
cally repeated) potential well of suitable shape. Suppose for a moment we are not 
adopting BvK boundary conditions, having thus only a genuinely isolated potential 
well. In this case the eigenstates can belong to two different classes: bound (local- 
ized) states, and scattering (delocalized) states. The distinction is a qualitatively 



clearcut one, and can be stated in several ways. One of them is to consider the 
second cumulant moment, or spread: 
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which is finite for bound states and divergent (when using appropriate normaliza- 
tions) for scattering ones. But if we study the same potential well within BvK, the 
qualitative distinction is lost: all states appear in a sense as "delocalized" since all 
wavefunctions ip(x) are periodic over the BvK period. And in fact the integrals in 
Eq. (8) become ill defined. 

The main issues therefore are: How do we distinguish between localized and 
delocalized states within BvK? In case of a localized state, how we actually measure 
the amount of localization? In the literature, such issues have been previously 
addressed by means of the participation ratio [15]. The complex number 3, whose 
phase provides the center of the distribution, Eq. (7), is our key to addressing 
localization: it is enough to consider its modulus. It has already been observed 
that |3| is bounded between and 1, and that I3I equals zero for an extremely 
delocalized state with |-^(a;)| 2 — 1/L. If we take instead an extremely localized 
state, with n\ oc (x) = 5(x), it is straightforward to get I3I = 1. It is therefore 
natural to measure localization by means of the negative of the logarithm of |3|: 
it is a nonnegative number, equal to zero in the case of extreme localization, and 
divergent in the case of extreme derealization. A glance at Eq. (6) yields: 



hence a natural expression for measuring the actual spread within BvK is: 



Having in mind again the eigenstates of a potential well, we can study the expression 
in Eq. (10) as a function of L. For a localized state, the shape of of ni oc (^) can be 
taken as L-independent for large L, hence Eq. (10) goes to a finite limit, which is 
the "natural" spread of the distribution. Quite on the contrary, for a delocalized 
state the distribution is smeared all over the (0, L) segment, preserving the norm 
over one period: therefore 3 goes to zero and the spread diverges in the large-L 
limit. 






(10) 



DIPOLE AND LOCALIZATION FOR MANY 

ELECTRONS 



So much about the one-electron problem: we are now going to consider a finite 
density of electrons in the periodic box. To start with, irrelevant spin variables will 



be neglected, and a system of spinless electrons in one dimension is considered. Even 
for a system of independent electrons, our approach takes a simple and compact 
form if a many-body formulation is adopted; BvK imposes periodicity in each 
electronic variable separately. Our interest is in studying a bulk system: N electrons 
in a segment of length L, where eventually the thermodynamic limit is taken: 
L — > oo, N — > oo, and N/L = n constant. 

We start defining the one-dimensional analogue of R of Eq. (1), namely, the 
multiplicative operator X = J2iLi x ii an d the complex number 

iN = (*\e i %*\*). (11) 

It is obvious that the operator X is ill-defined in our Hilbert space, while its 
complex exponential appearing in Eq. (11) is well defined. The main result of 
Ref. [4] is that the ground-state expectation value of the position operator is given 
by the analogue of Eq. (7), namely: 

(X) = |-Im In 3JV , (12) 

a quantity defined modulo L as above. 

The right-hand side of Eq. (12) is not simply the expectation value of an operator: 
it is the phase of it, converted into length units by the factor L/(2n). This phase 
can be called a single-point Berry phase, for reasons explained elsewhere [16,10,11]. 
Furthermore, the main ingredient of Eq. (11) is the expectation value of the multi- 
plicative operator e*T? x : it is important to realize that this is a genuine many-body 
operator. In general, one defines an operator to be one-body whenever it is the 
sum of N identical operators, acting on each electronic coordinate separately: for 
instance, the X operator is such. In order to express the expectation value of 
a one-body operator the full many-body wavefunction is not needed: knowledge 
of the one-body reduced density matrix p is enough: I stress that, instead, the 
expectation value of e l ^ x over a correlated wavefunction cannot be expressed in 
terms of p, and knowledge of the iV-electron wavefunction is explicitly needed. In 
the special case of a single-determinant, the iV-particle wavefunction is uniquely 
determined by the one-body reduced density matrix p (which is the projector over 
the set of the occupied single-particle orbitals): therefore the expectation value 
(X), Eq. (12), is uniquely determined by p. But this is peculiar to uncorrelated 
wavef unctions only. 

The expectation value (X) is extensive, as the dipole in Eq. (1). For the corre- 
sponding intensive quantity we borrow from Ref. [12] a useful notation: 

(x) c = (X)/N = -L-lm In 3jv , (13) 



where the subscript means "cumulant". The quantity (x) c goes to a well defined 
termodynamic limit, which is in fact proportional to the macroscopic polarization 



of the system. This result is proved in Ref. [9]; its three dimensional generalization 
is discussed in Refs. [10-12]. 

We stress that nowhere have we assumed crystalline periodicity. Therefore our 
definition of (x) c is very general: it applies to any condensed system, either ordered 
or disordered, either independent-electron or correlated. In the special case of a 
crystalline system, either interacting or noninteracting, the present approach can 
be shown equivalent to the previous formulations of polarization theory [1-4,10-12]. 

We are now ready to discuss electron localization in a condensed system: the 
present view is the one of Refs. [9,11], recently reexamined by Souza et al. [12], 
who also discuss its relationship to Kohn's localization [8]. This view is based 
on the modulus of 3 at, in full analogy with the previous Section about the single 
electron. We define therefore an intensive quantity, the second cumulant moment, 
by analogy with Eqs. (10) and (13): 

(- 2 )c = -^(^) 2 log|3v| 2 , (14) 

where again the notation is borrowed from Ref. [12]. This second moment is a very 
meaningful measure of electron localization in the electronic ground wavefunction, 
and enjoys two important properties: (1) When applied to a crystalline system 
of independent electrons, we recover an important gauge-invariant quantity which 
controls the Marzari-Vanderbilt localization [17]; (2) Even for more general sys- 
tems, correlated and/or disordered, (x 2 ) c assumes a finite value in insulators, and 
diverges in metals. Indeed in a metal the modulus of $n goes to zero in such a way 
that its phase is ill defined, and hence macroscopic polarization is ill defined as well. 
We have emphasized throughout this paper that one of the main phenomenological 
features differentiating insulators from metals is that the former materials sustain 
a nontrivial bulk polarization, while the latter do not. The complex number 3^ 
provides the key formal link between polarization and localization, via its phase 
and its modulus. 

NONINTERACTING ELECTRONS 

For independent electrons, we may write the many-body wavefunction $ as a 
Slater determinant of Bloch orbitals, but in the case of a metal not all the Bloch 
vectors in the reciprocal cell correspond to occupied orbitals: this fact is of over- 
whelming importance. We consider the simple case of one band in one dimension, 
whose Bloch vectors are illustrated in Fig. 2, imposing BvK boundary conditions 
over M cristal cells: L = Ma, where a is the lattice constant. We restore elec- 
tron spin here: we get an insulator if the number of electrons N equals 2M (filled 
band, top sketch), and a metal if N — M (half-filled band, bottom sketch). Both 
in the insulating and in the metallic case the iV-electron wavefunction is a Slater 
determinant of size N, built of N/2 doubly occupied spatial orbitals. Following the 




FIGURE 2. Discrete q vectors in the reciprocal cell in one dimension, where a BvK periodicity 
of M = 14 crystal cells has been chosen. Black circles correspond to occupied q's, and empty 
circles to unoccupied ones. Top: insulator, with N = 2M (filled band). Bottom: metal, with 
N = M (half-filled band). 



same algebra as in Refs. [4,10,11], the complex number $ N can be written in any 
case (insulator or metal) by means of the determinant of a matrix 

3w = (det S) 2 , (15) 

whose elements are 

and these elements are nonzero whenever s — s' + 1 . 

In the insulating case, owing to complete filling, both s and s' run over all the M 
values: for any given s, there is always one (and only one) s' such that the matrix 
element in Eq. (16) is nonzero. That means that in any row of the S matrix — 
whose size is M x M — there is one, and only one, nonvanishing element. Under 
these circumstances, the determinant factors as a product of M numbers: 

M-l 

det <S = H - / dx r qs+1 (^ q3 (x)^ x , (17) 
s =o a J 

where the identity ip qM (x) = ip qo { x ) is understood (periodic gauge). All the factors 
are nonvanishing, and the logarithm of 3^ is therefore a finite number. It can be 
shown that the N — > 00 limit of (x 2 ) c coincides with the spread of the optimally 
localized Wannier functions, as defined by Marzari and Vanderbilt [17]. 

The metallic case is very different. Since not all the q s vectors are occupied, the 
indices s and s' run over a subset of the M values (Fig. 2): the matrix S is of size 
M/2xM/2. There is one of the two q s at the Fermi level, for which the integrals in 
Eq. (16) are all vanishing, for any occupied s'. Therefore the matrix S has a row of 
zeros, and its determinant vanishes: its phase, and hence macroscopic polarization, 
is ill defined. The logarithm of 3^ is formally —00, and the spread (x 2 ) c diverges 
to +00. This is what we expected for a metal; the nontrivial fact is that it diverges 



(16) 



even at finite N, while on general grounds we only expected it to diverge in the 
thermodynamic (N — > oo) limit. 

So far, the compact and elegant expression of Eq. (14) has been proved to be 
appropriate to discriminate between insulators and metals only for a crystalline 
systems of independent electrons. For the — much more interesting — general case 
of a correlated and/or disordered system, we postulate that Eq. (14) performs the 
same task. The postulate is based on the general argument — much stressed above — 
about macroscopic polarization: well defined in insulators, ill defined in metals (as 
a bulk property). The correctness of this postulate has been verified by Resta and 
Sorella [9] for a one-dimensional model of a correlated crystal. Work on a model 
disordered system is in progress. Other very interesting discussions about the 
physical meaning of {x 2 ) c and its relationships to the insulating/metallic character 
of the system can be found in a paper of Souza et al. [12]. 



LOCALIZATION IN A MODEL CORRELATED SYSTEM 



We review here the very recent work of Resta and Sorella [9], where Eq. (14) 
is implemented for a one-dimensional two-band Hubbard model at half filling, 
intended to mimic an insulator having a mixed ionic/covalent character, and whose 
ground wavefunction is explicitly correlated. The macroscopic polarization of this 
model system was studied in Refs. [18,19]. To the present purpose, it is enough to 
study the centrosymmetric geometry: polarization is zero, the wavefunction is real, 
and the phase of 3 n is either or ir. 

The model has a very interesting behavior as a function of U: at small U it is a 
band insulator, while at a critical U c undergoes a transition to a Mott-like insulating 
phase. In the centrosymmetric geometry, $ N is a real number, which changes sign 
at U c . Its phase 7 {i.e. the single-point Berry phase [16,10,11]) jumps therefore by 
7r: it turns out that this occurrence is the main fact signalling the transition: the 
topological quantum number 7 /n can be used as an order parameter to identify the 
two different phases of the system [20,21]. At low values of U our model is a band 
insulator, while above U c is a Mott-like insulator. What happens to the second 
cumulant moment (squared localization length) (x 2 ) c as a function of U7 

The results are shown in Fig. 3 in terms of the dimensionless quantity 



V N = -JVlog |3jv 



(18) 



such that the squared localization length, Eq. (14) is; 




(27rn ) 2 n^oo 



N- 



(19) 



At U — the system is noninteracting, and the squared localization length (x 2 ) c 
coincides with the second moment of the (optimally localized) Wannier function of 
the occupied band. 




FIGURE 3. Dimensionless localization parameter, Eq. (18), as a function of U, after Ref. [9], 
where the effective value of TV is 800. The divergence of the localization length at the Mott 
transition is perspicuous. 

In the correlated case at U ^ no Wannier analysis can be performed; yet {x 2 ) c 
mantains its role as a meaningful measure of the localization of the electronic wave- 
function as a whole. The localization length increases with U below the transition, 
diverges at the transition point U c , and becomes localized again in the highly corre- 
lated regime, where in fact (x 2 ) c decreases with increasing U. Since the localization 
length remains finite at all values of U different from U c , the system is always in- 
sulating except at the transition point; the two insulating phases are topologically 
different and correspond to a qualitatively different organization of the electrons in 
the wavefunction [19]. 

What about the transition point? According to the previously stated viewpoint, 
the delocalized behavior implies a metallic character of the many-electron system. 
Notice that (x 2 ) c is a pure ground state property and apparently carries no in- 
formation about the excitation spectrum of the system. Yet we have explicitly 
verified — by exploiting the metastability of the Lanczos algorithm — that at the 
critical U value there is indeed a level crossing. At the transition value U c the 
ground state is twice degenerate and the lowest lying excitation (at constant N) 
has vanishing energy. 
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